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Magnetic reconnection during collisionless, stressed, X-point collapse was studied using kinetic, 
2.5D, fully electromagnetic, relativistic Particle-in-Cell numerical code. Two cases of weakly and 
strongly stressed X-point collapse were considered. Here descriptors weakly and strongly refer to 
20 % and 124 % unidirectional spatial compression of the X-point, respectively. In the weakly 
stressed case, the reconnection rate, defined as the out-of-plane electric field in the X-point (the 
magnetic null) normalised by the product of external magnetic field and Alfven speeds, peaks at 
0.11, with its average over 1.25 Alfven times being 0.04. During the peak of the reconnection, 
electron infiow into the current sheet is mostly concentrated along the separatrices until they defiect 
from the current sheet on the scale of electron skin depth, with the electron outflow speeds being 
of the order of the external Alfven speed. Ion infiow starts to defiect from the current sheet on 
the ion skin depth scale with the outfiow speeds about four times smaller than that of electrons. 
Electron energy distribution in the current sheet, at the high-energy end of the spectrum, shows a 
power law distribution with the index varying in time, attaining a maximal value of —4.1 at the final 
simulation time step (1.25 Alfven times). In the strongly stressed case, magnetic reconnection peak 
occurs 3.4 times faster and is more efficient. The peak reconnection rate now attains value 2.5, with 
the average reconnection rate over 1.25 Alfven times being 0.5. Plasma inflow into the current sheet 
is perpendicular to it, with the electron outflow seeds reaching 1.4 Alfven external Mach number 
and ions again being about four times slower than electrons. The power law energy spectrum for the 
electrons in the current sheet attains now a steeper index of —5.5, a value close to the ones observed 
near X-type region in the Earth's magneto-tail. Within about one Alfven time, 2% and 20% of the 
initial magnetic energy is converted into heat and accelerated particle energy in the case of weak and 
strong stress, respectively. In the both cases, during the peak of the reconnection, the quadruple 
out-of-plane magnetic field is generated, hinting possibly to the Hall regime of the reconnection. 
These results strongly suggest the importance of the collisionless, stressed X-point collapse as an 
efficient mechanism of converting magnetic energy into heat and super-thermal particle energy. 

PACS numbers: 52.35.Vd; 96.60.Iv; 52.65.Rr; 45.50.Dd; 96.60.pf; 96.60.qe 



I. MOTIVATION OF THE STUDY 

Magnetic reconnection is an important physical pro- 
cess, which serves as one of the possible ways of convert- 
ing energy stored in the magnetic field into heat and non- 
thermal, accelerated, motion of plasma particles. This 
process operates virtually in all extra-galactic, stellar, 
solar, space and laboratory plasmas with varied degree 
of importance. For example, in solar and stellar flares 
magnetic reconnection plays a key role. In addition, it 
can be one of the main contributing factors to solar coro- 
nal heating problem amongst other mechanisms such as 
wave dissipation. As far as plasma heating is concerned, 
in the coUisional regime, be it Tokamak plasma or solar 
corona (which perhaps is better described as a collision- 
less medium), Spitzer resistivity is oc Te where Te 
is the electron temperature. Thus, in becoming hotter, 
plasma essentially starts to behave as a superconductor, 
i.e. further heating (increase in temperature) is impeded 
due to decrease of the resistive properties. In this con- 
text, a direct conversion of magnetic energy into heat 
via reconnection seems rather attractive. As far as the 
particle acceleration is concerned, 50-80 % of the energy 
released during solar flares is converted into the energy of 
accelerated particles. Thus, knowledge of the details of 



particle acceleration via magnetic reconnection (which is 
deemed as a major mechanism operating solar and stellar 
flares) is also important. 

The main aspects of the reconnection process can be 
classified as whether it is resistive (collisional) or colli- 
sionless; spontaneous or forced; and/or steady or time- 
dependent. Each of these aspects have been studied ex- 
tensively with varied level of progress. For example, both 
steady and time-dependent, resistive reconnection pro- 
cesses are very well studied [see e.g. Ref. [l| and references 
therein] , while it is only recently some progress has been 
made in the study of collisionless reconnection [see e.g. 
Ref. 0, Q and references therein] . This can be explained 
by the fact that resistive reconnection is simpler in its na- 
ture. Changing structure (connectivity) of magnetic field 
lines, which essentially is the reconnection, requires some 
mechanism (dissipation) to break the frozen-in condition. 
In the case of resistive reconnection it is r/j term. In fact, 
rate at which reconnection proceeds is given by an inflow 
velocity into the dissipation region Vin ~ {5 / IS)va, where 
5 and A are width and length of the region, and va is 
Alfven speed. In the simplest resistive, steady model of 
Sweet-Parker 5 ~ rj^^"^ and A ~ L, where L is the macro- 
scopic system size. In the collisionless regime however, 
other terms in the Ohm's law, such as electron inertia 
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term, Hall term and electron pressure tensor become far 
more important than the usual r]j term on which resistive 
reconnection models rely. Each of these terms has an as- 
sociated spatial scale with them, i.e. scale at which they 
start to play a dominant role. For example, electron iner- 
tia term is associated with the electron skin depth c/ujpe, 
Hall term is with the ion skin depth c/ujpi^ while electron 
pressure tensor is with ion Larmor radius [see e.g. p. 87 
in Ref.d or p. 42 in Ref.llij]. 

One of the main outcomes of recent research in colli- 
sionless reconnection can be summarised by GEM [1, Q 
and Newton [6| reconnection challenges. These consid- 
ered the same physical system: Harris-type equilibrium 
with anti-parallel magnetic field, relevant to geomagnetic 
tail application, with finite initial magnetic perturbations 
in the case of GEM challenge and time-transient inho- 
mogeneous, driven inflow of magnetic flux in the case 
of Newton challenge. The novelty of the approach was 
to use different numerical codes [MHD, Hall MHD, Hy- 
brid and Particle-in-Cell (PIC)], in order to pin-point es- 
sential physical mechanism that facilitates the reconnec- 
tion. These works established that as long as dispersive 
whistler waves are included (these can only appear if one 
allows for the different dynamics for electrons and ions), 
the rate at which reconnection proceeds does not change, 
irrespective of which term breaks the frozen-in condition. 

Yet another interesting analytical result corroborated 
by the numerical simulations is that in 2D steady re- 
connection, Petschek reconnection rate functional form 
remains the same when the Hall term is included (in the 
Hall MHD model), but what changes is the length from 
the X-point to the start of slow mode shocks Q. In 
MHD, this length is the half width of the resistive re- 
gion, Ln, but when Hall term is included the length is 
replaced by Lji + c/cupi j7|. This important result essen- 
tially provides an analytical scaling for the coUisionless 
reconnection (when rj is zero) for such physical system. 
Thus, future kinetic studies that use PIG simulation need 
to corroborate it. Ref.Q already tested this scahng law 
via MHD and Hall MHD simulation. 

Importance of the electron inertia in the X-point with 
finite out-of-plane guide magnetic field has been inves- 
tigated in Ref. Particularly noteworthy result was 
that when the normalised coUisionless electron skin depth 
[c/{ujpeL)] exceeds the dimensionless resistive length scale 
(5"^/^), where S is the Lindquist number, the energy in 
a shear Alfven wave approaching an X-point is rapidly 
transformed into plasma kinetic energy and heat. As- 
suming solar coronal parameters, the normalised electron 
skin depth is 10~^ (assuming L = 10 Mm and c/ujpe = 0.1 
m), while dimensionless resistive length scale is 3 x 10~^ 
(in the corona typically 5* = 10^^). Thus, the proposed 
in Q may well be effective. 

Previous works on particle acceleration mostly focused 
on test-particle type calculations of the particle trajec- 
tories in different 2D [§, [l^, IHl 113 or more recently 
3D [IB, [3 magnetic reconnection configurations. 

In such approach feedback on reconnection electromag- 



netic (EM) fields from motion (spatial redistribution) of 
charged particles is ignored. Our approach does not suf- 
fer from this drawback, as we use fully electromagnetic, 
relativistic PIC numerical code in which EM fields are 
calculated at each step from the spatial distribution of 
the charges. 

According to Ref.fl], resistive time-dependent recon- 
nection other than well known tearing mode can be split 
into two main classes: X-type collapse, which was first 
considered by Ref.fl^l (a decade before tearing mode 
was discovered) and Petschek-type theory developed by 
Ref.[ll]. Previous work on X-point collapse is well de- 
scribed in chapter 7.1 of Ref.[l|]. Also see more recent 
work on the subject [l^- A good example of combina- 
tion of analytical and numerical work on magnetic recon- 
nection at stressed X-type neutral points can be found 
in Ref.jl^]. Boundary conditions used in Ref.jl^l, where 
such that they did not allow for flux or mass flow through 
the boundary. The studies with closed boundaries are 
physically justified by being isolated systems, whereas 
some of those with open boundaries lead to misleading 
results: for example, in an open system a potential X- 
point can collapse due to inflow of energy from outside. 

Initial analytical work on this topic considered un- 
bounded self-similar solutions. These indicated that 
-^2(0,0,*), the out-of-plane electric held at the magnetic 
null, which is the measure of reconnection rate, tends 
to infinity as t 00 [e.g. Fig. (7.3) from Ref.[l|]. The 
main outcomes of the previous, stressed X-point collapse 
in the resistive MHD (in the case of low-resistivity and 
low-beta) can be summarised as following: 

(i) The X-point collapse is different depending whether 
initial stressing is weak or strong. In the weak case 
the average reconnection rate scale as 1/ ln(?7) , where 
f] — 1/S is the dimensionless resistivity [77 = •ql{VAoL) = 
1/5], while strong case it is independent of fj. 

(ii) There is an issue related to efficiency of the process: 
Ref. 21] showed that for fast reconnection to occur (3 < 
j^o.565^ otherwise pressure in the current sheet chokes off 
the collapse process. For solar coronal conditions, fj ~ 
10^^^ this requires f3 < 10~^, which is too low to match 
the observed values circa 0.01 — 0.001. However, this 
shortcoming seems to be alleviated by inclusion of non- 
linear effects, i.e. case of strong perturbations (strong 
stressing) [2l|. If the compression is sufficiently large, 
magnetic pressure of the imploding wave expels trapped 
gas in the current sheet from its ends, allowing it to thin 
and faster reconnection to occur. 

Motivation of the present study is five-fold: 

(i) Naturally, different boundary and initial conditions 
produce different scaling of reconnection rate e.g. with 
resistivity (in the case of resistive reconnection). To our 
knowledge X-type collapse by a uniform stress (compres- 
sion) along one of the axis, as described in chapter 2.1 
m 

Ref.51, has not been investigated numerically neither 
in the case of resistive (MHD) reconnection, nor in the 
coUisionless regime; 

(ii) This type of stress (compression) is likely to occur 
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FIG. 1: Magnetic field fine geometry, (a) A portion of a 
standard Solar flare model with added compressions on sides. 
The bottom part mimics the footpoints of a coronal loop, 
with an X-point on top (inside the dashed box). Two arrows 
indicate compression direction. A dashed box is the region 
which our study tries to mimic by uniform stress in one direc- 
tion [see (c)]. (b) Magnetic field configuration for stress pa- 
rameter Q = 1, which is stable (no X-point collapse occurs), 
(c) Magnetic field configuration for stress parameter a > 1. 
Such X-point collapses under high magnetic pressure/tension 
forces. 



the sides the considered magnetic configuration is sta- 
ble [Fig. [TJb)]. Interestingly, Ref. [231] investigated what 
happens when a fast magnetosonic shock wave associ- 
ated with a coronal mass ejection collides obliquely with 
a coronal streamer with a stable current sheet. Their set 
up is somewhat analogous to ours, but much more vio- 
lent, as we only consider X-point compression, while they 
blasted it with a Alfven Mach 6 shock. 

The initial magnetic field configurations used are 

B.^^y, By^^a'x, =0, (1) 

where _Bo is magnetic field intensity at the distance L 
from the X-point for a = 1.0, L is the global external 
length-scale of reconnection, and a is the stress parame- 
ter [see chapter 2.1 in Ref. [l|]. In addition, the uniform 
current is imposed at t = in the ^-direction, 

j, = ^(a^-l). (2) 

This current appears because of the stress imposed on 
magnetic field lines. 



in the framework of a flare model (see details below) . 

(iii) To our knowledge the most of previous collision- 
less reconnection studies considered anti-parallel, Harris 
type magnetic field configuration, which is more relevant 
to Earth magneto-tail application. For solar and stellar 
fiares X-type configuration is more relevant. 

(iv) Some of the models of coronal mass ejections [l^l 
which use motion of photospheric footpoints as a driver, 
end up with situation physically similar to stressed X- 
point collapse (strictly speaking Y-points occur there). 

(v) We aimed to investigate properties of the acceler- 
ated particles in the current using self-consistent electro- 
magnetic fields, as a further extension of the test-particle 
type approach. 



II. SIMULATION MODEL 

A. Stressed X-point reconnection model 

Fig. [H represents magnetic field line geometry consid- 
ered in this paper. In order to study stressed magnetic 
reconnection, we consider a magnetic X-point collapse 
which may naturally occur during solar flares [see sketch 
in Fig.[T]Ja)]. Such stressed X-point may occur e.g. if pho- 
tospheric footpoints of the coronal loops move towards 
each other (e.g. pushed by convective motions) or some 
compression from the sides in the corona forces them to 
do so. We study dynamics of such stressed X-point by 
means of kinetic, 2.5D, fully electromagnetic, relativistic 
Particle-in-Cell numerical code. We focus attention on 
the local region inside the dashed box which is is the re- 
gion that our study tries to mimic by uniform stress in 
one direction [see Fig. [BJc)]. If there is no stress from 



B. PIC simulation code 

The simulation code used here is 2.5D relativistic and 
fully electromagnetic PIC code, modified from the 3D 
TRISTAN code fH]. In this code, both the electron 
and ion dynamics are described as particles. Equation 
of motion for each particle is solved using self consistent 
electromagnetic fields. 
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dt eo 



dB 

a7 = -^^^' 



(3) 



(4) 



(5) 



(6) 



where i?, -B, ja, Vsi and rgi are electric and magnetic 
fields, current density, particle velocity and position, re- 
spectively. The subscript s represents species of plasma, 
that is s = e for electron and s — i for ion. The subscript 
i indicates i-th particle index. The other quantities, g^, 
TOs, c, eo, are charge and mass of a plasma particle, speed 
of light, vacuum permittivity, respectively. In addition 
to above Eqs. ©-(H]), V • i? = and V • pe/eo 
must be satisfied initially (pe is charge density, which is 
taken as 0). The latter two conditions are automatically 



4 



A - 




h .... I .... I .... I .... I .... I ... . 
5 10 15 20 25 30 



FIG. 2: Time evolution of the reconnected magnetic flux dif- 
ference. The horizontal axis is time normalised by the ion 
cyclotron frequency, uod- The solid curve is the magnetic 
flux difference, A-i/i, between the O and X lines, which is 
normalised by Boc/uJpi, using our PIC code. Open squares 
correspond to published GEM challenge values. 



satisfied at all times due to the nature of the numerical 
scheme used 24|. Electromagnetic fields which satisfy 
Eqs. (O-® update particle velocity through the equa- 
tion of motion Eq. ([3]), which in turn updates particle 
positions Eq. Repeating this procedure many times 
mimics plasma particle dynamics in self-consistent elec- 
tromagnetic fields. 

We started by testing our code, reproducing published 
GEM challenge results [4;,^]. Fig.[2]shows time evolution 
of the reconnected magnetic flux difference between and 
O and X-lines. Open squares correspond to the published 
GEM challenge values. We gather from this graph that 
match between our simulation results and that of GEM 
challenge [1, Q is good. We could also reproduce other 
figures from Refs.jJQ, thus we were reasonably confident 
to start a new X-point collapse simulation, which is the 
subject of the present paper. 

The length of the system in two dimensions is = 
Ly = 400A (this is excluding, so called, ghost cells needed 
for updating the boundary conditions), where A = 1 is 
the simulation grid size corresponding to electron Debye 
length, Xo — vte/i-^pe — lA (vte is electron thermal ve- 
locity and LOpe is electron plasma frequency) . The global 
external length-scale of reconnection is set L = 200 A. 
Number density is fixed at no = 100 electron-ion pairs 
per cell. Both electrons and ions are uniformly dis- 
tributed throughout the system, hence total number is 
1.6 million pairs. 

Zero-gradient boundary conditions are imposed both 
on the electric and magnetic fields in x- and y- 
directions. Also, tangential component of electric field 
was forced to zero, while normal component of mag- 
netic field was kept constant, both at the bound- 
ary. This ensures that there is no magnetic flux 



through the simulation boundary, i.e. the system is iso- 
lated and — J^j^ By{x' , — L, t)dx' — /^^ By{x' , L, t)dx' + 

j':^B,{~L,y\t)dy' + j":^ B,{L,y' ,t)dy' = 0, at every 
time step. The latter sum is the magnetic flux on the 
boundary. Reflecting boundary conditions are imposed 
on particles in both the x- and y-directions. The latter 
ensures there is no mass flow across the boundary. 

The simulation time step is LUpgAt = 0.05. Ion to elec- 
tron mass ratio is rrii/me = 100. Electron thermal ve- 
locity to speed of light ratio is vte/c — 0.1. Electron 
and ion skin depths are c/tOpe = lOA and c/ujpi = lOOA, 
respectively. Electron cyclotron frequency to plasma fre- 
quency ratio is ujce/^^pe — 1.0 for magnetic field intensity, 
B = Bq. The latter ratio is indeed close to unity in the 
solar corona, while it is much bigger than unity in the 
Earth magnetosphere. Electron and ion Larmor radii 
are vtej^ce = lA and vu/oJci = lOA, where va is the ion 
thermal velocity. The temperatures of ions and electrons 
are the same, Te = Ti. At the boundary the plasma 
13 = 0.02 and Alfven velocity, Vao/c = 0.1. Naturally 
these vary across the simulation box as the background 
magnetic field is a function of x and y. 

In what follows, for the visualisation purposes, all spa- 
tial coordinates will be normalised by electron skin depth 
c/cjpe, while time is normalised by the inverse of plasma 
electron frequency ui^^- 



III. SIMULATION RESULTS 

We investigated three X-point collapse cases with dif- 
ferent stress parameters, a = 1.00 (stable), 1.20 (weakly 
stressed) and 2.24 (strongly stressed). For a — 1.00, we 
confirmed that the system is stable at least for t < 500, 
and no magnetic reconnection takes place. Ref.[2i| has 
shown that when perturbations (exerted stress on an X- 
point) are small e ^ {1 — a^) < fj then the (average) re- 
connection rate scales as l/ln(^), while if they are large 
e — (1— a^) > fj then the reconnection rate is inde- 
pendent of fj. This difference in behaviour points to a 
different relative importance of the physical mechanisms 
in action. The two different (weakly and strongly com- 
pressed) cases studied in this paper attempt the same 
approach as in Ref.[2l[. However, one should realise that 
we use vastly different physical description. Ref . [2l| uses 
resistive MHD, while our PIC (fully kinetic) approach 
is collisionless, and, in turn, in our case ^ = 0. How- 
ever, it has been observed before, that often scattering 
of plasma particles from the magnetic field lines, plays 
effective role of collisions. E.g. Refs.[25l. [26j have shown 
that the Alfven wave dissipation in the collisionless case 
(using the same PIC, kinetic code) follows scaling law 
which is the same as in resistive MHD case [23|- At the 
same time, this is not to say that we are confident that 
the effective particle scattering off the magnetic fields, 
which mimics resistive (coUisional) effects, is the mecha- 
nism that breaks down the frozen-in condition here. The 
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FIG. 3; (a) Time evolution of the reconnection rate, defined 
as the out-of-plane electric field in the X-point, Ez{0, 0, t) nor- 
malised by the product of external magnetic field and Alfven 
speeds (both at the boundary for a = 1.00) for a = 1.20 
(solid line) and a = 1.00 (dotted line), (b) Time evolution of 
total current density, jz, in the X-point for a — 1.20. Note 
that final simulation time t — 250 corresponds to I/ta ~ 1.25, 
with the latter being the Alfven time. 



issue of which term in the generalised Ohm's law is re- 
sponsible for the reconneetion in the coUisionless stressed 
X-point collapse will be studied separately. 



A. Weakly stressed X-point case (a = 1.20) 

It should be noted that a, the measure of stress, is 
essentially the aspect ratio of the compression, i.e. since 
the limiting magnetic field lines are given hy y = ax, 
a — 1.2 means 1.2-1.0=0.2 i.e. 20% compression of the X- 
point in x-direction. Similarly, a = 2.24 case represents 
124% compression, i.e. 1.24 times stressed. 

Also, because change in a means change of strength 
of the magnetic field on the boundary, the external 
Alfven speed also changes. E.g. we fixed Alfven veloc- 
ity, Vao/c = 0.1 for a = 1, but for a > 1, Vaq/c = 
V(l + a^)/2 X 0.1. This for a = 1.20, Vao/c = 0.124, 
but for a = 2.24, Vao/c = 0.36. 

Fig- EI a) shows time evolution of the reconnection rate, 
defined as the out-of-plane electric field in the X-point, 
Sz(0,0,t) normalised by the product of external mag- 
netic field and Alfven speeds (both at the boundary), i.e. 
by Eo = VaqBq, for a = 1.20 (sohd fine) and a = 1.00 
(dotted line). Because usually PIC simulation suffers 
for large (thermal) noise especially towards high energy 
ends of the particle distribution function, a boxcar aver- 
age scheme with width 150 mesh points was applied for 
smoothing the line-data. Mind that this does not alter 
absolute value of the peak aX t — 170. Our normalisa- 
tion of the electric field is such that effectively Fig. [31^ a) 
shows the Alfven Mach number of the inflowing plasma. 



Ma, i.e. E/Eq = v/Vao = Ma- This follows from a typ- 
ical estimate E ~ vB^ ■^i MaVaoBq [e.g. Ref.|10|]]. We 
gather from this figure that for the case of a = 1 no i?2 
is generated, as such configuration is stable and no X- 
point collapse or reconnection occurs (dotted line). On 
contrary, for the case of a = 1.2 we see that the recon- 
nection rate peaks at about 0.11 at t = 170 which is the 
same as I/ta = 0.85. Here ta is the Alfven time - the 
distance from the outer boundary to the X-point (20) di- 
vided by the Alfven speed at the boundary (V^o = 0.1c). 
We also calculated the average reconnection rate based 
on Fig. [21 a). It is defined as 

tf Jo Eo 

where tf = 250 is the final simulation time. This yields 
Eav =0.04, a relatively small value. 

Fig.[3ljb) shows time evolution of total current density, 
jz, in the X-point for a = 1.20 (solid line) and a = 1.0 
(dotted line). The total current density is normalised by 
the initial value, jo = noevdo, Vdo is the drift velocity. 
As expected a = 1 case stays current free throughout 
the simulation, while a = 1.2 produces a peak current of 
jz/jo = 15 and then subsequently decays off. 

It should be noted that we have performed one addi- 
tional numerical run with a = 1.2, but without imposing 
initial jo current prescribed by Eq.(2). This is because 
e.g. one might expect that in a X-point above the arcade 
of loops in the solar corona, only compression (stress) of 
the magnetic field from the two sides is likely. Such per- 
turbation violates the equations at t — 0, this results in 
a large spike of magnitude 0.22 at f = 5 in the equivalent 
version of Fig. [3Ka) (not presented here). Otherwise, for 
t > 5 there was no noticeable difference between the cases 
with and without imposing the initial current. Thus, in 
what follows we only discuss cases with the current im- 
posed a.t t — 0. 

Fig. 2] shows time evolution of spatial distribution of 
total current density, j^, in the x-y plane at (a) t — 0, (b) 
100, (c) 170 and (d) 250 for a = 1.20. The current flows 
uniformly in the z-direction Fig. IH^a) , with drift velocity 
is Vdo/c — 0.022. We assumed that half of the drift ve- 
locity comes from electrons and another half from ions, 
i.e. Vdeo = O.Swdo and Vdto = O-^Vdo- Fig. g^b) shows 
that the current begins to focus in the X-point because 
of high magnetic pressure/tension. The current peaks in 
the current sheet attaining jz/jo = 15 at i = 170, as seen 
in Fig.m^c). In the later phases, the current sheet tends 
to decay [Fig. (Hd)]. This flgure is useful for visualis- 
ing the spatial dimensions of the current sheet. In turn, 
this enables to make the following useful estimate: From 
Fig. d^c) (peak current time snapshot) we gather that 
the width of the current sheet (in horizontal direction) 
is about electron skin-depth 5 ~ 1.4, while its length (in 
vertical direction) is about A ~ 10. To be precise, we 
actually measured 6 and A by more accurate means: we 
looked at hue plots of jz{x,0) and jz{0,y) respectively, 
at t = 170, and measured appropriate half-width of the 
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FIG. 4: Time evolution of the spatial distribution of total 
current density, jz, in the X-Y plane at (a) t = 0, (b) 100, 
(c) 170 and (d) 250 for a — 1.20. The total current density is 
normalised by the initial value, jo = noevdo- 



FIG. 5: Time evolution of the spatial distribution of the out- 
of-plane magnetic field, Bz, at (a) t — 0, (b) 100, (c) 170 
and (d) 250 for a = 1.20. The magnetic field intensity is 
normalised by the initial value, Bq. 



jzix,0) and jz{0,y) peak (in both a: and y-directions) . 
The ratio of the two gives the inflow Alfven Mach number 
Ma — S/ A ~ 0.14. This is very close to the peak value of 
E/Eo = v/Vao = Ma^ 0.11 [from Fig.^a)]. We thus 
measured the peak reconncction rate by two independent 
means directly from the simulation [using Fig. [31(a)] and 
using the steady reconnection formula Ma — 5/ A [e.g. 
Eq.(3.1) from Ref.^Si]]. The close match points to the fact 
that nearly steady reconnection regime is achieved in the 
vicinity of the peak aX t = 170. 

Fig. [5] shows the time evolution of spatial distribution 
of the out-of-plane magnetic field, B^, at (a) t = 0, (b) 
100, (c) 170 and (d) 250 for a = 1.20. Time of each panel 
in Fig. [5] corresponds to those of Fig. SI respectively. As 
can be seen in Fig. [DJa) initially there is no out-of-plane 
magnetic field present. Then, Fig. [5I^b) shows that the 
quadruple magnetic field structure appears, which is the 
most pronounced &t t = 170 [Fig. EI^c)] . This structure 
is a well known signature of a magnetic reconnection in 
the Hall MHD regime. Ref. ^28] showed that as long as 
one allows for different electron and ion dynamics (two- 
fluid description) i.e. when Hall term is non-zero such 
quadruple out-of-plane magnetic field is generated. Bz 
attains a value of \Bz\/Bq = 0.18 at t = 170 and then 
subsequently decays, Fig. O^d), (as well as the current 
sheet). 

In Fig. [6] we show dynamics of individual magnetic field 
lines. We tried to trace dynamics of several magnetic field 
lines in order to visualise the reconnection process. We 
can clearly see that magnetic field lines come towards 
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FIG. 6; Time evolution of the magnetic field lines in the x-y 
plane at (a) t = 0, (b) 100, (c) 170 and (d) 250 for a = 1.20. 
The several lines are plotted with different intensity, \B\/ Bq = 
1.55, 1.60, 1.65 and 1.70 on the boundaries. The left and right 
field lines in the simulation box are distinguished with solid 
and dotted line styles, respectively, to visualise magnetic field 
reconnection clearly. 
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FIG. 7: Electron (a) and ion (b) flow pattern around X-point o.oooi o.ooio o.oioo o.iooo i.oooo 

at t = 170 for a = 1.20. Note that here Vao/c = 0.124. E/(m.„c^) 



each other in x-direction, reconnect, and move apart in 
y-direction. 

Fig. [7] shows electron (a) and ion (b) flows at the peak 
time of the reconnection. It is rather instructive to see 
that this figure in effect corroborates the sketch from 
Ref.[l|, [see their Fig. (3.1)]. In particular it shows that 
the electron and ion flow are clearly separated on the 
two different spatial scales - electron skin-depth and ion 
skin-depth, note that since here the mass ratio is 100, 
c/uopi = lOc/cjpe)- The noticeable difference is caused be- 
cause they considered initial background magnetic fleld of 
Harris- type (anti-parallel), while we study X-point. Nat- 
urally, in our case electron inflow into the current sheet 
is mostly concentrated along the separatrices until they 
deflect from the current sheet on the scale of electron skin 
depth, with the electron outflow speeds being of the order 
of the external Alfven speed 0.13c. Ion inflow starts to 
deflect from the current sheet on the ion skin depth scale 
with the outflow speeds about four times smaller (0.03c) 
than that of electrons. As argued by Ref.[2§] it is dif- 
ference in these two flows which generates the observed 
quadruple out-of-plane magnetic fleld. 

Fig.[8]shows the local electron energy distribution func- 
tion near the current sheet at t = (dashed curve) and 
t = 250 (solid curve) for a = 1.20. We performed numer- 
ical fit to the high energy part of the distribution function 
{E > O.OStoeoc^ = 41 keV) and found that the electron 
distribution has a power law form i.e. in particular the 
best fit is provided by / = dN/dE oc E~'^-^ (straight 
solid line). In general, the direct observational evidence 
for the relation between reconnection and particle ac- 
celeration is hard to obtain because acceleration length 
scale is of the order of ion skin-depth. However, recently 
it became possible to directly measure electron energy 
spectrum in the vicinity of X-type region in the Earth's 
magneto-tail [13] ■ They found that from about 2 keV to 
about 200 keV power law index was between -4.8 and -5.3 
(changing in the course of time of the observation). It is 
interesting to note that both our power law index (-4.1) 
and the energy (0.2 x meoc^ ~ 0.2 x 511 ~ 100 keV) at 
the high end of the spectrum are close to the observed 



FIG. 8: The local electron energy spectrum (distribution 
function) near the current sheet at f = (dashed curve) 
and t = 250 (solid curve) for a — 1.20. Data is pro- 
duced using region near the maximum current sheet gener- 
ation (—2 < X < 2, —8 < y < 8). The vertical axis shows 
the number of electrons with a particular energy. The lat- 
ter is shown on the horizontal axis (normalised by electron 
rest energy, nieoc^)- Straight solid line is our best fit, while 
dash-dotted and dash-triple-dotted lines are taken from the 
observations (see text). 

values by Ref. [1^ . Note that since the parameters of our 
simulation are for the solar corona {ujce/^pe — 1-0), and 
are not even totally realistic for that (e.g. to simplify 
numerical simulation our mi/ni^^ Vte/c^ and Vao/c ra- 
tios are not entirely realistic), the absolute values of our 
simulation quoted in keV should be taken with caution 
(when comparing them to Earth's magneto-tail results). 
Yet we are confident in the correctness of the obtained 
power law indexes. Also, noteworthy fact is that the 
whistler wave turbulence, based on Fokker-Planck equa- 
tion for the electron distribution function, subject to a 
zero-flux boundary condition (same as ours by a chance) , 
is producing [3l| similar power law energy spectrum. Re- 
call that it is standing whistler waves are thought to be 
mediating reconnection in the Hall regime. 

B. Strongly stressed case (a = 2.24) 

In this subsection we present results of numerical sim- 
ulation for the case of a = 2.24, which is regarded as a 
case of strongly stressed X-point. 

Fig. [HI shows the numerical simulation results as in 
Fig. [3] but for a — 2.24. To avoid repetition we omit- 
ted a = 1 case [dotted line in Fig. [3] . We gather from 
panel (a) that now reconnection rate attains value of 2.5 
at i = 45 (0.225t^). In the strongly stressed case, in 
differ the weakly stressed one, Ez rebounds and an oscil- 
lation is established. We calculated average reconnection 
rate based on Fig.[51[a) using definition from Eq.(7). The 
result is Eav = 0.5, a value ~ 10 times bigger than in the 
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FIG. 9: The same as in Fig. [3] but for a = 2.24. 



weakly stressed case, indicative of more efficient recon- 
nection. Also, much stronger current density, jz/ ja = 55, 
is generated compared to the weak case [panel (b)] and 
peak of the reconnection occurs 170/45 = 3.8 times ear- 
lier. Similar oscillations of current and electric field were 
used as a mechanism for interpreting the peculiar hard 
x-ray (> 25 keV) solar flare, which is believed to be pro- 
duced by a non-thermal electron beam [s^, [ssj . 

Fig. [TO] panels (a) and (b) show spatial distribution of 
the current at the peak of the reconnection t = 45, and 
after the current sheet decayed, t = 70, [the first bounce - 
see Fig.|9]. The noticeable difference from panels (c) and 
(d) in Fig.|3]is that now current sheet becomes longer and 
thinner [the same conclusion reached in Ref.[21|, in that 
large perturbations (strong compression) yields current 
sheet thinning and onset of more efficient reconnection]. 
As in the weakly stressed case, we also did the following 
calculation: Using data from Fig. [TOTa) (peak current 
time snapshot) the width and length of current sheet was 
estimated by looking at half-width of the jz{x,0) and 
jz{0,y) peak (in both x and y-directions). The result is 
width, S ~ 0.4, and length, A ~ 17 The ratio of the two 
gives the inflow Mach number Ma ~ J/A ~ 0.02. This is 
125 times smaller that peak reconnection rate of 2.5 from 
Fig. [Ifa) and 25 time smaller than average reconnection 
rate of 0.5 from the same figure. This discrepancy can 
only be attributed to the fact that formula A/ ~ (5/A only 
applies to the steady reconnection process. In the case 
of a — 1.2 the reconnection process was relatively steady 
thus the latter formula proved a good match, but for 
a = 2.24 the process is too dynamic and thus Ma — S/A 
as a measure of reconnection fails. Fig. [TO] panels (c) and 
(d) again show familiar quadruple out-of-plane magnetic 
field structure, but now it is much more elongated. The 
panels (e) and (f) where we trace individual magnetic 
field lines at two different times provide another proof 
that the reconnection takes places and that the current 
sheet is now much longer and thinner. 

In Fig. [TT] is an analog of Fig. [7] but for a = 2.24. 
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FIG. 10: Spatial distribution of total current density, jz, at 
(a) t = 45 and (b) t = 70 for a = 2.24. The total current 
density is normalised by the initial value, jo = noevdo- The 
same for out-of-plane magnetic field, B^, which is normalised 
to Bo at (c) t ^ 45 and (d) t = 70. Panels (e) and (f) 
show individual magnetic field lines with different intensity, 
\B\/Bo = 2.0, 3.5, 4.0 and 4.5 on the boundaries at t = 45 
and (d) t = 70 (f) respectively. The left and right field lines 
in the simulation box are distinguished with solid and dotted 
line styles, respectively. 



Noteworthy difference from the weak case is that the 
plasma inflow into the current sheet is perpendicular to 
it, with the electron outflow speeds reaching external 
Alfven Mach number of 0.5/0.36 =1.4 [note arrow length 
in Fig. Illf a)]. and ions again being about four times 
slower than electrons. Mind that now Vao/c = 0.36. 

Fig. [12] shows the local electron energy distribution 
function in the current sheet at t = (dashed curve) 
and t = 250 (sohd curve) for a = 2.24. Note that the 
dashed curves (corresponding to t — 0) in Fig. [5] and 
Fig.[T2]are different. This is due to the fact that different 
a's, mean different initial jo current [see Eq.(2)], which 
contribute to the calculation of the distribution function. 
We performed numerical fit to the high energy part of the 
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FIG. 11: The same as in Fig.Qbut for a = 2.24. Here t = 45. 
Note that here Vao/c — 0.36. 
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FIG. 12: The same as in Fig.|8]but for a = 2.24. Here t = 250. 
Data is produced using a region near the maximum current 
generation ( — 1 < a; < 1, —16 <y< 16). 

distribution function {E > 2.4meoc^ = 1.226 MeV) and 
found that the electron distribution has the power law 
index of -5.5 which is quite close to the observed power 
law range of -4.8 and -5.3 [13] ■ Note that now attained 
energies are much higher 4 x Weoc^ — 0.2 x 511 keV ~ 2 
MeV (at the high end of the spectrum). Again, we note 
that values quoted in keV and MeV should be taken with 
caution [see discussion of Fig. [5] above] . 



C. Energetics of the reconnection process 

We now try to estimate the efficiently of the X-point 
collapse by looking at the energetics of the process. 
In particular in Fig. [T3] we plot the magnetic energy, 
Esit) ^ J J{B^2 + By2 + B^2)/{2fio)dxdy,foi two cases 
a — 1.2 (solid line) and a — 2.24 (dashed line). The nor- 
malisation in each case is Eb{0)- Note that the latter is 
different in the both cases because different stress means 
different initial magnetic field. We gather from this graph 
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FIG. 13: The magnetic energy, EB{t) = / f{Ba:2 + By2 + 
Bz2)/{2fio)dxdy, normalised on its initial value, Eb{0), for 
two cases a = 1.2 (solid line) and a = 2.24 (dashed line) 
versus time. Note that ujp^t = 250 corresponds to 1.25 Alfven 
times). 



that in the case of weak stress, 2% of the initial magnetic 
energy (which is a dominant part of the total energy be- 
cause our plasma beta is 0.02) is converted into heat and 
energy of super-thermal (accelerated) particles. In the 
case of large perturbations (strong compression) 20 % in 
the initial magnetic energy is converted into heating and 
particle acceleration. 



IV. CONCLUSION AND DISCUSSIONS 

We studied magnetic reconnection during coUisionless, 
stressed, X-point collapse using kinetic, 2.5D, fully elec- 
tromagnetic, relativistic Particle-in-Cell numerical code. 
We investigated two cases of weakly and strongly stressed 
X-point. Here by weakly and strongly we mean 20 % and 
124 % unidirectional spatial compression of the X-point, 
respectively. The reconnection rate, defined as the out- 
of-plane electric field in the magnetic null normalised by 
the product of external magnetic field and Alfven speeds, 
peaks at 0.11 (at 0.85 Alfven times), with its average 
over 1.25 Alfven times being 0.04. We found that dur- 
ing the peak of the reconnection, electron inflow into the 
current sheet is mostly concentrated along the separatri- 
ces until they deflect from the current sheet on the scale 
of electron skin depth, with the electron outflow speeds 
being of the order of the external Alfven speed. Ion in- 
flow starts to deflect from the current sheet on the ion 
skin depth scale with the outflow speeds about four times 
smaller than that of electrons. Electron energy distribu- 
tion in the current sheet, at the high energy end of the 
spectrum, shows a power law distribution with the index 
which varies in time, attaining a maximal value of -4.1 at 
the final simulation time step (this corresponds to 1.25 
Alfven times). 
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The obtained results in the strongly stressed case show 
that the magnetic reconnection peak occurs about 3.4 
times faster (at 0.25 Alfven times) and is more efficient. 
The peak reconnection rate now attains value 2.5 (at 0.25 
Alfven times), with the average reconnection rate over 
1.25 Alfven times being 0.5. Plasma inflow into the cur- 
rent sheet is perpendicular to it, with the electron out- 
flow seeds reaching 1.4 Alfven external Mach number and 
ions again being about four times slower than electrons. 
The power law energy spectrum for the electrons in the 
current sheet attains now a steeper index of -5.5. This 
is close to the typical observed power law indexes in the 
vicinity of X-type region in the Earth's magneto-tail [13] ■ 

The reconnection rate versus time figures in both cases 
indicate that the reconnection has bursty, time-transient 
behaviour. In the considered two cases, 2% (in weakly 
stressed) and 20% (in strongly stressed) of the initial 
magnetic energy is converted into heat and energy of ac- 
celerated particles, respectively, both within about one 
Alfven time. This is somewhat similar to the previous 
resistive MHD analog [2l| of the present study, in that 
small perturbation (weak stress) initiates X-point col- 
lapse, but then reconnection process is choked off by the 
gas pressure. At present, it is unclear, however, whether 
we seem to observe similar behaviour in our simulation. 
After all, e.g. Petschek mechanism reconnection can be 
choked off by a different (other than pressure) mechanism 
[when the diffusion-region inflow magnetic field gets too 
smafl [!]. 

We also found that in the both cases, during the peak 
of the reconnection, the quadruple out-of-plane magnetic 
field is generated. This most is likely to suggest that Hall 
regime of the reconnection takes place [28j . 

In addition to the fundamental interest of converting 
magnetic energy into heat and the energy of accelerated 
particles, these results are significant for e.g. solar coro- 



nal heating problem. It has been estimated that even 
if only 2% of the magnetic energy in the solar corona is 
converted into heat, the coronal heating problem would 
be solved. Also, it is known that the typical resistive dif- 
fusion time in the corona is 10^^ Alfven times (10^ years), 
while the time-transient phenomena such as flares (one 
of the possible candidates of coronal heating) occur on 
time scales of 10 - 100 Alfven times. Our results for the 
strongly stressed case suggest that 20 % of initial mag- 
netic energy can be released in just 1 Alfven time. On one 
hand, one has to realise, that the obtained results are for 
an X-point configuration, and naturally, the bulk of solar 
corona is not made of solely of X-points. On the other 
hand, it is also known that heating of the active regions 
would provide circa 82 % of the coronal heating budget 
[34]. In turn, given the fact that X-points are common 
occurrence in the active regions, our results seem to be 
of importance for solving the coronal heating problem. 

It important to stress that the obtained main result 
of the 20% conversion of the initial magnetic energy into 
heat and energy of super-thermal particles, within about 
one Alfven time, is obtained in the coUisionless regime, 
thus this result does not suffer from any uncertainty in 
the anomalous resistivity as in the case of resistive MHD 
reconnection. 



Acknowledgments 

Authors acknowledge use of computational facilities 
available through UKMHD consortium lead by Univer- 
sity of St. Andrews. This research was funded by 
the United Kingdom's Science and Technology Facilities 
Council. 



[1] E. Priest and T. Forbes, Magnetic reconnection: MHD 
theory and applications (Cambridge University Press, 
2000). 

[2] D. Biskamp, Magnetic reconnectton m Plasmas (Cam- 
bridge University Press, 2005). 

[3] J. Birn and E. R. Priest, Reconnection of magnetic fields: 
magnetohydrodynamics and coUisionless theory and ob- 
servations (Cambridge: Cambridge University Press, 
2007). 

[4] J. Birn, J. F. Drake, M. A. Shay, and et al., J. Geopliys. 

Res. 106, 3715 (2001). 
[5] P. L. Pritchett, J. Geophys. Res. 106, 3783 (2001). 
[6] J. Birn, K. Galsgaard, M. Hesse, and et al., Geophys. 

Res. Lett. 32, L06105 (2005). 
[7] T. D. Arber and M. Haynes, Physics of Plasmas 13, 2105 

(2006). 

[8] K. G. McClements, N. Shah, and A. Thyagaraja, J. 

Plasma Phys. 72, 571 (2006). 
[9] G. E. Vekstein and P. K. Browning, Physics of Plasmas 

4, 2261 (1997). 



[10] P. K. Browning and G. E. Vekstein, J. Geophys. Res. 

106, 18677 (2001). 
[11] V. V. Zharkova and M. Gordovskyy, Astrophys. J. 604, 

884 (2004). 

[12] B. Hamilton, L. Fletcher, K. G. McClements, and et al., 

Astrophys. J. 625, 496 (2005). 
[13] S. Dalla and P. K. Browning, Astron. Astrophys. 436, 

1103 (2005). 

[14] S. Dalla and P. K. Browning, Astrophys. J. Lett. 640, 
L99 (2006). 

[15] P. Browning and S. Dalla, Mem. Soc. Astron. Italiana 
78, 255 (2007). 

[16] Y. E. Litvinenko, Astron. Astrophys. 452, 1069 (2006). 

[17] J. W. Dungey, Phil. Mag. 44, 725 (1953). 

[18] V. S. Semenov, M. F. Heyn, and I. V. Kubyshkin, Sov. 

Astron. 27, 660 (1983). 
[19] G. Vekstein and N. Bian, Astrophys. J. 632, L151 

(2005). 

[20] L. Ofman, P. J. Morrison, and R. S. Steinolfson, Astro- 
phys. J. 417, 748 (1993). 



11 



[21] A. N. McClymont and I. J. D. Craig, Astrophys. J. 466, 

487 (1996). 

[22] T. G. Forbes and E. R. Priest, Astrophys. J. 446, 377 
(1995). 

[23] J. 1. Sakai and Y. Tanaka, Astron. Astrophys. 468, 1075 
(2007). 

[24] O. Buneman, Computer Space Plasma Physics, Simula- 
tion Techniques and Software ed. H. Matsumoto & Y. 
Ornura (Terra Scientific, 1993). 

[25] D. Tsiklauri, J. 1. Sakai, and S. Saito, New J. Phys. 7, 
79 (2005). 

[26] D. Tsiklauri, J. 1. Sakai, and S. Saito, Astron. Astrophys. 

435, 1105 (2005). 
[27] J. Heyvaerts and E. R. Priest, Astron. Astrophys. 117, 

220 (1983). 



[28] D. A. Uzdcnsky and R. M. Kulsrud, Physics of Plasmas 

13, 2305 (2006). 
[29] M. A. Shay, J. F. Drake, B. N. Rogers, and et al., J. 

Geophys. Res. 106, 3759 (2001). 
[30] M. 0ieroset, R. P. Lin, T. D. Phan, and et al., Phys. 

Rev. Lett. 89, 195001 (2002). 
[31] C.-Y. Ma and D. Summers, Geophys. Res. Lett. 25, 

4099 (1998). 

[32] V. M. Nakariakov, C. FouUon, E. Verwichtc, and ct al., 

Astron. Astrophys. 452, 343 (2006). 
[33] L. Ofman and L. Sui, Astrophys. J. 644, L149 (2006). 
[34] M. J. Aschwanden, A. Winebarger, D. Tsiklauri, and 

et al, Astrophys. J. 659, 1673 (2007). 



